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We compute the bispectrum for the IRAS PSCz catalog and find that the galaxy distribution dis- 
plays the characteristic signature of gravity. Assuming Gaussian initial conditions, we obtain galaxy 
biasing parameters 1/bi — 1.20^o'jg an d 62/6? — —0.42 ±0.19, with no sign of scale-dependent bias 
for k < 0.3 h/Mpc. These results impose stringent constraints on non-Gaussian initial conditions. 
For dimensional scaling models with x% statistics, we find N > 49, which implies a constraint on 
primordial skewness B3 < 0.35. 
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Introduction. — What is the nature of the distribution 
of matter in the universe? Does the galaxy distribu- 
tion follow the mass distribution on large scales? Our 
ability to constrain cosmological parameters from galaxy 
clustering is only as good as our understanding of bias, 
that is, the relation between the galaxy distribution and 
the underlying dark matter distribution. In the cur- 
rent structure formation paradigm, most of the matter 
in the universe is dark, interacting only through grav- 
ity, and cold. Small fluctuations, generated in the early 
universe, grow gravitationally once the universe becomes 
matter-dominated, leading to the formation of large-scale 
structure. Non-linear gravitational instability (GI) the- 
ory makes specific predictions for the statistical proper- 
ties of the mass distribution, e.g., how non-Gaussian fea- 
tures arise from Gaussian initial conditions. The extent 
to which the galaxy distribution follows these predictions 
can be used to put constraints upon galaxy bias (and 
thus upon the physics of galaxy formation) and upon the 
degree of non-Gaussianity of the primordial fluctuations 
(which constrains the physics of the generation of fluctu- 
ations in the early universe). 

In this Letter, we report and interpret measurements of 
the bispectrum for the IRAS PSCz galaxy redshift survey 
[jl| . Details of the method and earlier results can be found 
in two previous papers. In the first B, we developed 
a framework to obtain reliable constraints upon galaxy 
biasing and primordial non-Gaussianity that takes into 
account the effects of survey geometry and sampling. In 
the second |J , we computed the bispectrum of the galaxy 
distribution of earlier IRAS redshift catalogs. Due to 
its larger volume and higher sampling density, the PSCz 
survey yields the best evidence supporting GI and the 
strongest constraints to date upon galaxy biasing and 
primordial non-Gaussianity. 

The Bispectrum. — In linear perturbation theory (PT), 
the statistical properties of the primordial fluctuations, 
Gaussian or non-Gaussian, are preserved. If Gaussian, 



they can be fully characterized in Fourier space by their 
power spectrum P(k), 

(S m (ki)S m (k 2 )) = (2vr) 3 ( 5 D (fc 1 + fe 2 ) P m (ki) , (1) 

where S m (k) is the Fourier transform of the density con- 
trast Sp/p. Assuming a linear and deterministic galaxy 
bias, 6(x) = bi5 m (x), measurement of the galaxy power 
spectrum in redshift surveys in the linear regime gives 0] 



P(k,p) = (1 + 2/V + /?V) (M 8 ) 2 P!(fc), 



(2) 



where (3 w il^/bi, p is the cosine of the angle between 
the wavector and the line of sight, as is the rms density 
fluctuation amplitude in spheres of radius 8 h' 1 Mpc, and 
Pi (k) is the linear mass power spectrum scaled to erg = 1- 
Equation (0) illustrates the degeneracies between the pa- 
rameters ft m , 61, and as' measurement of the quadrupole 
to monopole ratio of P(k,p) via redshift distortions, for 
example, yields a measurement of /?, whereas the over- 
all amplitude of the power spectrum gives in addition an 
estimate of bias- Equation (J2J) involves the small angle 
approximation, however for the parameter combinations 
measured here, the results will be unchanged even if the 
approximation is dropped. 

Additional information is needed to break this degen- 
eracy between f2 m , b±, and as] for example, non- linear 
effects in pairwise velocities J5|,^| and in density-velocity 
comparisons can be used to constrain different com- 
binations of f2 m , bi, and as- We use here measurements 
of the bispectrum (see also [|]-[l2| ) . The nonlinearities in 
the equations of motion of GI induce a characteristic pat- 
tern of three-point and higher order correlations Jp^JllI . 
For Gaussian initial conditions, the induced bispectrum 
-B123 = B(ki,k 2 ,k 3 ), 



(S m (ki)5 m (k 2 )S m (k 3 )) 



(27r) 3 Mfci 



fc 2 + k 3 ) B123, 
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scales as B(k) oc [P(k)] 2 \ it is thus useful to examine the 
reduced bispectrum amplitude Q, 
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To leading order in perturbation theory for Gaussian ini- 
tial conditions, to an excellent approximation GI predicts 
that Q is independent of time and thus of eg and f2 m ; Q is 
a function of triangle shape with a behavior that depends 
only on the spectral index, n c s(k) — d\nPi(k)/dlnk. 

We observe galaxies, not mass. At scales larger than 
the characteristic scale of galaxy formation, it is reason- 
able to assume that the galaxy density is given by a local 
transformation of the underlying mass density field; that 
is, for the smoothed density S = b\8 m + 562^ + • • • lo- 
cally, in which case the galaxy bispectrum amplitude (Q) 
is related to the matter amplitude (Q m ) [PJ§ as 



Or, 

h 



(5) 



This shows how the galaxy bispectrum can be used to 
break the degeneracies present in linear theory: since 
Q m is independent of <r 8 and Q m , measuring Q for dif- 
ferent triangles and comparing it to the predicted Q m 
based upon the measured spectral index gives separate 
values for the bias parameters 61 and !)2 |. In redshift 
space (where distances are measured by recession veloci- 
ties using the Hubble law) there are small corrections to 
Eq. (|) that can be computed; for b\ = 1, 62 = —0.5, 
and f2 m = 0.3, underestimates the exact result in the 
large-scale limit by AQ ks 0.14, independent of config- 
uration p^ |. Thus, it does not affect b\, and given the 
observational uncertainties, it only minimally alters the 
inferred 62- 

To estimate the galaxy power spectrum and bispec- 
trum, we follow Feldman, Kaiser, and Peacock |ll| in 
using the optimal weighted transform of the difference 
between the data and a tailored synthetic catalog 
In addition we construct mock catalogs using second- 
order Lagrangian PT (2LPT), which agrees with TV-body 
simulations at the scales we consider g]. The speed of 
2LPT allows us to construct a large number of mock 
catalogs which take into account the selection function, 
survey geometry, and redshift-space mapping. This is 
essential to estimate accurately the mean, errors, covari- 
ance, and the distribution of Q for different triangles, 
required to compute the non-Gaussian likelihood func- 
tion The 2LPT mock catalog initial conditions are 
based on a ACDM model with Q m = 0.3, fl A = 0.7 (A 
stands for the cosmological constant); the linear power 
spectrum shape parameter T — Cl m h — 0.21 is consis- 
tent with PSCz power spectrum measurements p7| , ^8[ , 
and the adopted normalization, erg = 0.7, corresponds to 
ct| — 0.84 in redshift space, consistent with the observed 
normalization for IRAS galaxies fl9|]. 
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FIG. 1. Q from the PSCz survey vs. k = ki for triangles 
with ratio of sides fe/fci separated by angle 9. Panels show 
bands in triangle shape, from nearly collinear (upper left) to 
open (lower right). Triangles, squares, and circles show results 
for ki/ki = 0.25-0.5, 0.5-0.75, 0.75-1.0. The dashed lines in 
each panel show the median and 68% range of the points. 

Data and Analysis. — The IRAS PSCz survey (!]] is 
based upon the IRAS Point Source Catalog with addi- 
tions to achieve the best possible uniformity over the sky. 
Out of the 15411 IRAS galaxies across 84% of the sky, 
we take the 13180 galaxies in the range 20 /i -1 Mpc < 
R < 500 h' 1 Mpc, with galactic latitude \b\ > 10°, and 
IRAS 60 micron flux / 60 > 0.6 Jy. Our results appear 
graphically in Figures 1 and 2. The bispectrum is defined 
for closed triangles in Fourier space. We characterize a 
given triangle by its ordered sides k\ > ki > k%. In Fig. 1 
we show Q from PSCz as a function of k\ for triangles 
with two sides separated by angle cos# = (fei • ^2), for 
four ranges of 0. Since the width of the survey window 
function is « 0.01 h Mpc -1 , we confine our quantitative 
analysis to triangle sides k% > 0.05 h Mpc -1 , for which 
deconvolution using the narrow window approximation 
is accurate [||. The median of Q decreases with increas- 
ing 9, as predicted by GI Jl3| (see Fig. 2). 

Figure 2 shows Q from PSCz as a function of angle 
for triangles with 0.2 < k\ < 0.4 h Mpc -1 and with two 
sides of ratio k^jkx = 0.4-0.6. The solid curve shows the 
prediction for Q in redshift space, obtained from averag- 
ing many 2LPT mock catalogs as detailed above. Sym- 
bols show results for bands in k\, as identified in the 
caption. The dependence of Q on triangle shape pre- 
dicted by GI is clearly seen in the data. We employ a 
likelihood analysis as in to quantify the constraints 
that this places on galaxy biasing and primordial non- 
Gaussianity. For Gaussian initial conditions, the dashed 
curve shows the 2LPT ACDM predictions with the best- 
fit values I/61 = 1.20, b 2 /b\ = -0.42, as given in Table 1. 
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FIG. 2. Q from the PSCz survey vs. # for triangles 
with 0.2 < ki < 0.4 h Mpc -1 and with sides of ratio 
fe/fci = 0.4 — 0.6 separated by angle 8. The solid curve shows 
Q in redshift space averaged over many 2LPT realizations of 
the ACDM model. The dashed curve shows the prediction 
for ACDM with the best-fit bias parameters 1/bi — 1.20, 
b 2 /b\ = -0.42. 

Figure 3 shows the likelihood contours and marginal- 
ized distributions for l/b\ and &2/&1 for the PSCz sur- 
vey, assuming Gaussian initial conditions and including 
triangles with k < /c max = 0.3 ft. Mpc -1 , where the upper 
limit is imposed by both nonlinearity and discreteness. 
We have varied the maximum wavenumber allowed in 
the analysis, fc max = 0.2, 0.25, and 0.3 ft. Mpc -1 , (see Ta- 
ble 1.) The maximum likelihood values are essentially 
independent of k max , though the error bars decrease as 
/: max is increased. There is no evidence that the bias pa- 
rameters depend on k. These results for b\ and 62 are 
insensitive to the values of erg and Q m adopted in the 
mock catalogs [||||, but they do depend upon the as- 
sumed value for T, since it determines the predicted Q m 

Si- " 

If galaxy biasing is local in Lagrangian, rather than Eu- 
lerian space (imposed on the initial, rather than final con- 
ditions), then Q in Eq.(|^) has an additional term which 
depends upon configuration pfj| . We have repeated the 
likelihood analysis for such a model and found that the 
likelihood ratio for the Eulerian to Lagrangian maximum 
likelihood is £e/£-l = 1-6. This implies that the be- 
havior of the bispectrum for PSCz galaxies is better de- 
scribed by an Eulerian bias scheme. 

TABLE I. Likelihood Estimates of Bias Parameters 



parameter 


^max — 0.2 


^max — 0.25 


&max — 0.3 


I/61 
62/6? 


1 1 c+0.36 
*— La -0.37 

-0.38±°;£ 


1 1 1 +0.22 
1 - 1L -0.2S 

-0 29+ ' 27 


1 20+°- 18 
-0 42+ - 19 




FIG. 3. Likelihood contours for bias parameters b\ and 
62 from the PSCz survey for fc ma x = 0.3 h Mpc -1 , assuming 
Gaussian initial conditions and a linear power spectrum with 
F = 0.21. The contours denote 68%, 90% and 95% confidence 
levels. 

The PSCz bispectrum provides strong constraints 
upon non-Gaussian initial conditions. As an extension 
of the results in || , we consider \n statistics as a gen- 
eral example of dimensional scaling models, where the 
linear PT n-point correlation function £„ oc ^ 2 . For 
N = 1, this corresponds to the predictions of some infla- 
tionary models with isocurvature perturbations pl| , p2| ; 
as N — > 00 the model becomes effectively Gaussian, 
and for a fixed power spectrum (taken to fit that of 
PSCz) the primordial bispectrum obeys Q 1 oc A^ -1 / 2 
p3[ . In addition, there are two other contributions that 
partially cancel that scale as A^ -1 j23j. The amplitude 
and scaling of the PSCz bispectrum constrains A" to be 
considerably larger than unity. Thus, we use Q(N) ~ 
Qg + (Q x 2 ~ Qg)^^ 1 / 2 , and similarly for the errors, 
where Qg and Q x 2 are obtained from mock catalogs for 
Gaussian and N = 1 initial conditions, respectively. In 
addition, we assume that Q-eigenmodes and eigenvalues 
are well approximated by the Gaussian ones in the large 
N regime. We construct a likelihood for A?" by marginal- 
izing over the bias; from this we obtain N > 49 (22) at 
68% (95%) CL. Since the primordial dimensionless skew- 
ness is B3 = 2.46 for a \\ field [M i we conclude that the 
PSCz bispectrum constrains B3 < 0.35(0.52). 

Discussion. — Our findings are consistent with previ- 
ous studies from sparser and smaller IRAS surveys ||, 
but using the PSCz catalog reduces the statistical and 
systematic errors and allows us to place the strongest 
constraints to date upon nonlinear bias and primordial 
non-Gaussianity. Our main results are the following: 
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(1) The dependence of the bispectrum upon configu- 
ration shape is the clearest direct evidence to date that 
non-linear mode couplings are as specified in GI. 

(2) The linear bias parameter we find, 1/bx = 



1.20±^|, agrees well with b x = 0.84 ±0.28 [f24| obtained 
by a joint analysis of the PSCz power spectrum, cos- 
mic microwave background (CMB) anisotropy, and clus- 
ter abundances. It is also consistent with that inferred 
for APM galaxies |^5| using the shape of the two-point 
correlation function, &apm = 1-11 ± 0.22 [fZq , taking 
into account the relative bias &apm/^pscz = 1-15 [ p7[ . 
It is in less good agreement but still consistent with 
frpscz = 1.10 ±0.23 (95% CL) found from a concordance 
between PSCz and recent CMB results [p8|. 

-0.19 
-0.18' 



(3) We obtain a nonlinear bias bijb\ = —0.48. 



which is not yet measured in any other way. A negative 
62 means that the number of IRAS galaxies does not rise 
as quickly as the mass density contrast; this is consistent 
with the observation that infrared-selected galaxies are 
increasingly under-represented in massive clusters. 

(4) That both the scale and shape dependence of the 
bispectrum are consistent with GI acting on Gaussian 
initial conditions limits the possibility of primordial non- 
Gaussianity. In particular, our study of \n initial condi- 
tions constrains N > 49 and the primordial skewness to 
be B 3 < 0.35. 

(5) Our results for the bispectrum are consistent with 
the galaxy skewness S3 measured for the same catalog 
p9[ and are similarly insensitive to the flux cut: we find 
identical results for / 60 > 0.8 Jy. Further, the PSCz data 
is incomplete beyond 420 h' 1 Mpc [Q, however, since we 
look at k > 0.05, this should not affect the bispectrum 
signal. 

(6) When coupled with measurements of the power 
spectrum redshift distortions, which determine f3 — 
fl° m 6 /bi ~ 0.4 ± 0.12 for the PSCz survey @2§], our 
measurement of b\ from the bispectrum allows one to 
break the degeneracy between linear bias and fl m : we 
find fi TO = 0.16 ± 0.1, indicating a low-density universe. 
We stress that, unlike other methods, this result is based 
only on the clustering properties of PSCz galaxies assum- 
ing Gaussian initial conditions, independent to a large 
extent of external determinations of erg . 

Due to its almost complete sky coverage and large vol- 
ume, the PSCz survey is the best survey available to 
study the bispectrum at large scales. This will soon 
change with the advent of larger redshift surveys such as 
2dF and SDSS, which will determine the bispectrum even 
more precisely and thus test the GI framework, nonlinear 
galaxy bias, and primordial non-Gaussianity to unprece- 
dented accuracy. 
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